NASA/TM— 2007-214690 



An Introduction to Thermodynamic Performance 
Analysis of Aircraft Gas Turbine Engine Cycles 
Using the Numerical Propulsion System Simulation 
Code 


Scott M. Jones 

Glenn Research Center, Cleveland, Ohio 


March 2007 


NASA STI Program ... in Profile 


Since its founding, NASA has been dedicated to the 
advancement of aeronautics and space science. The 
NASA Scientific and Technical Information (STI) 
program plays a key part in helping NASA maintain 
this important role. 

The NASA STI Program operates under the auspices 
of the Agency Chief Information Officer. It collects, 
organizes, provides for archiving, and disseminates 
NASA’s STI. The NASA STI program provides access 
to the NASA Aeronautics and Space Database and its 
public interface, the NASA Technical Reports Server, 
thus providing one of the largest collections of 
aeronautical and space science STI in the world. 
Results are published in both non-NASA channels and 
by NASA in the NASA STI Report Series, which 
includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant phase 
of research that present the results of NASA 
programs and include extensive data or theoretical 
analysis. Includes compilations of significant 
scientific and technical data and information 
deemed to be of continuing reference value. 
NASA counterpart of peer-reviewed formal 
professional papers but has less stringent 
limitations on manuscript length and extent of 
graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies that 
contain minimal annotation. Does not contain 
extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, often 
concerned with subjects having substantial 
public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific and 
technical material pertinent to NASA’s mission. 

Specialized services also include creating custom 

thesauri, building customized databases, organizing 

and publishing research results. 

For more information about the NASA STI 

program, see the following: 

• Access the NASA STI program home page at 
http://www.sti.nasa.gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA STI Help Desk 
at 301-621-0134 

• Telephone the NASA STI Help Desk at 
301-621-0390 

• Write to: 

NASA Center for AeroSpace Information (CASI) 
7115 Standard Drive 
Hanover, MD 21076-1320 


NASA/TM— 2007-214690 



An Introduction to Thermodynamic Performance 
Analysis of Aircraft Gas Turbine Engine Cycles 
Using the Numerical Propulsion System Simulation 
Code 


Scott M. Jones 

Glenn Research Center, Cleveland, Ohio 


National Aeronautics and 
Space Administration 


Glenn Research Center 
Cleveland, Ohio 44135 


March 2007 


This report is a formal draft or working 
paper, intended to solicit comments and 
ideas from a technical peer group. 


Level of Review. This material has been technically reviewed by technical management. 


Available from 


NASA Center for Aerospace Information 
7115 Standard Drive 
Hanover, MD 21076-1320 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22161 


Available electronically at http://gltrs.grc.nasa.gov 


An Introduction to Thermodynamic Performance Analysis of 
Aircraft Gas Turbine Engine Cycles Using the Numerical 
Propulsion System Simulation Code 

Scott M. Jones 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Abstract 

This document is intended as an introduction to the analysis of gas turbine engine cycles using the 
Numerical Propulsion System Simulation (NPSS) code. It is assumed that the analyst has a firm 
understanding of fluid flow, gas dynamics, thermodynamics, and turbomachinery theory. The purpose of 
this paper is to provide for the novice the information necessary to begin cycle analysis using NPSS. This 
paper and the annotated example serve as a starting point and by no means cover the entire range of 
information and experience necessary for engine performance simulation. NPSS syntax is presented but 
for a more detailed explanation of the code the user is referred to the NPSS User Guide and Reference 
document (ref. 1). 


The Gas Turbine Engine 

The Brayton (or Joule) cycle best describes the operation of an air-breathing gas turbine engine. The 
engine consists of three main components — a compressor, burner, and turbine. For aircraft propulsion the 
cycle is open, meaning the turbine exhaust is vented to the atmosphere rather than recirculated back 
through the compressor. Chemical energy provided by the fuel is converted to useful work in one of two 
ways: the turbine can produce additional mechanical (shaft) power; or, a nozzle can provide a momentum 
increase in the working fluid (air) to produce thrust. Due to the complex interactions between the engine 
components as well as large changes in environmental operating conditions, computer codes are required 
to predict the performance of all but the most simple, ideal engine. NPSS is one such code and its object- 
oriented nature enables nearly any conceivable engine architecture to be accurately modeled. 


Basic NPSS Input File Syntax 

The basic format of an NPSS input file is text-based. This text can be broken down into three main 
categories: creation of objects, assignment of values to variables, and commands (or function calls). In 
object-oriented codes the creation of a particular instance of an object type is called instantiation. The 
input file is read (or parsed) sequentially so normally an object will be created/instantiated then values 
assigned to its pertinent variables and the process repeated with the next object. Once all objects have 
been declared then commands or function calls are made to the code. Because NPSS allows user-defined 
functions it is important to distinguish between creation of a function and actually calling or executing 
that function. 

When assigning values to variables it is important to understand variable “scope” which, loosely 
defined, is the region where a variable exists. In order to assign a value to a variable while outside its 
scope, any objects above that scope will have their name prepended to the variable name with a “dot” 
separator (e.g., object_l.object_2.variable_name = 10). See the “Variable Scope” section of the NPSS 
User Guide and Reference for a more detailed explanation. Finally, remember that NPSS is case sensitive 
and all commands and variable assignments must end with the semicolon character. 
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Example Model Input Structure 

Input files for NPSS may have any reasonable name or extension; a simple turbojet model consisting 
of several files is listed in appendix B. Appendix A is a detailed explanation of this example model. 
Because input files are read sequentially, models typically have a layout reflecting steps a through g 
outlined below: 

a) Comments 

Because input files may have different names, it is good practice to place a comment box at the 
beginning of every file explaining what that file does or the data it contains. The double slash (//) 
comments everything to the end of that single line, and multiple lines may be commented by 
using the slash and asterisk (/* ... */) at each end of the comment. 

b) Specify Thermodynamic Package 

The thermodynamic package to be used is specified first. Many NPSS objects cannot be created 
unless this has been done, and this is one of the few times a command will occur near the 
beginning of an input file. The thermodynamic package determines fluid properties (viscosity, 
enthalpy, thermal conductivity, etc.) from known fluid state conditions (pressure and temperature, 
for example). The command to specify the thermodynamic package is 
setThermoPackage(“name”); where ‘name’ is either ‘GasTbT, ‘Janaf , ‘allFueT, or ‘CEA’. 
Janaf is well suited to account for mixtures of air and water vapor typically found in a gas turbine 
environment. CEA is an implementation of the CEA code (ref. 2) which can account for many 
chemical species and interactions but at significant cost in computational time and stability. It is 
recommended that ‘Janaf be used unless chemical equilibrium calculations are required or 
another package is deemed more appropriate. 

c) Physical Components and Configuration (Linkages) 

Next, the model is actually configured by instantiating objects which represent the engine 
physical components and linking them to form the model structure. These objects are Elements 
and Subelements such as inlets, compressors, ducts, and turbines. Elements can be instantiated in 
any order, but a layout reflecting the way flow travels through the engine is recommended. Also, 
each Element is given a unique name by the user so that, for example, a model may have more 
than one compressor. This “building block” approach allows nearly any practical engine 
architecture, from turboshafts to recuperated turbofans, to be accurately simulated as a set of 
interconnected components. 

Once all the Elements are declared, their outlet and inlet fluid ports are connected by the 
linkPorts command. This determines the actual layout of the engine (e.g., the compressor exit is 
connected to the burner entrance, the burner exit is connected to the turbine entrance, the turbine 
exit is connected to the nozzle entrance, etc.). Similarly, shaft ports are also connected by the 
linkPorts command to represent mechanical linkages in the engine. 

The format of the command to connect Element ports is linkPorts(“elementl.outlet_port”, 
“element2.inlet_port”, “link_name”); where Tink name’ is a unique name for the li nk . 

d) Solver Objects (Independents and Dependents) 

In order to attain a physically valid solution, an NPSS engine model must at a minimum obey 
conservation of mass, energy, and momentum; the NPSS solver ensures this when properly 
configured. Furthermore, the solver can also alter independent variable values (e.g., design fan 
pressure ratio) so that a desired condition is achieved (e.g., engine net thrust equals 10000 lb). 
Independents and Dependents are the most common solver objects used in NPSS. Some Elements 
in NPSS already have Independents and/or Dependents associated with them, and others can be 
created by the user. The creation of an Independent looks like this: 
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Independent myindepname { 

varName = “actualvariablename”; 

} 

where ‘my indep name’ is a unique user-specified name for the Independent (because there are 
likely to be many Independent objects in the model) and “actual variable name” (double quotes 
are necessary) is the actual name of the variable whose value will be altered when this 
Independent is active. Note the use of curly braces to set the scope of the variable ‘varName’; an 
alternative is 

Independent my indep name; 

myindepname.varName = “actualvariablename”; 

The creation of a Dependent is similar: 

Dependent mydepname { 

eqrhs = “variable_parameter_l”; 
eqlhs = “variable_parameter_2”; 

} 

where ‘my dep name’ is a unique user-specified name for the Dependent (because there are 
likely to be many Dependent objects in the model). The variables ‘eq_rhs’ and ‘eq_lhs’ can be 
thought of as the right and left hand sides of an equation, or a condition to be satisfied. The 
strings “variable_parameter 1” and “variable_parameter_2” (double quotes are necessary) can be 
variable names, numerical values, or even calculations. As examples: 

eqrhs = “enginethrust”; 
eqlhs = “10000”; 

eqrhs = “enginethrust”; 
eqlhs = “desiredthrust”; 

eqrhs = “enginethrust”; 

eq lhs = “desired_thrust*2.0 - 5000.0”; 

Note that in these examples ‘engine thrust’ and ‘desired thrust’ must be recognized NPSS 
variables (although they can be user-created variables). 

e) Set Up Solver (Make a Set of Independents and Dependents Active) 

An engine model may contain any number of Independents and Dependents, but they are not 
automatically added to the solver upon their creation. Independents and Dependents are added to 
the solver by the solver.addIndependent(“my_indep_name”) and 

solver.addDependent(“my_dep_name”) commands. Conversely, active solver objects may be 
made inactive by the solver.removeIndependent(“my_indep_name”) and 
solver.removeDependent(“my_dep_name”) commands. The remove command only makes 
solver objects inactive, it does not delete them; they can be made active again later in the model 
simulation. 

When the run() command is issued, the NPSS solver will execute the engine model to 
convergence if possible. This means that each Element and Subelement will perform its 
calculate() function; and, any and all active Independent variable values will be altered to satisfy 
all of the active Dependent conditions. Elements are likely to have their calculate() function 
executed many times as the model iterates to a solution. 
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Note that although any number of Independents and Dependents can be created, the number of 
active Independents must always equal the number of active Dependents for any particular case at 
the time of the run() command. Finally, the user may send to the screen a listing of the active 
Independents and Dependents at any point in a simulation by the following commands: 

cout « solver .list( “Independent”, TRUE ) « endl; 
cout « solver .list( “Dependent”, TRUE ) « endl; 

This can be most helpful when evaluating the behavior of a model simulation. 

f) Run the Case 

Once the solver has been set up to provide the desired analysis, the run() command is issued. 

Data may then be written to one or more files immediately after the engine analysis case is run. A 
single input file will almost always run many cases; for each separate case a few variable values 
are changed (such as altitude and flight Mach number) and steps “e” and “f ’ are repeated. 

g) Output Format Specification (DataViewers) 

By default, no output is generated by NPSS; the user must either use the cout command (which 
sends data to standard output) or create one or more DataViewer objects which send output to a 
file (or an output stream). A DataViewer uniquely determines which variables to output as well as 
the look and format of the output. Each DataViewer object in the model is an object with a unique 
name and therefore the user can output data to them whenever and only as desired. To output the 
data to a file the command dataviewer_name.display() is issued, where ‘dataviewername’ is 
the unique name of that DataViewer. Each DataViewer also has an output stream object 
associated with it; this output stream specifies the actual file where the data is written (the 
‘filename’ string variable). There are several types of DataViewer objects and their use is 
explained further in the NPSS User Guide and Reference. 

With the exception of DataViewers, which can be specified almost anywhere in the input file as long 
as it is prior to their display command being issued, the creation of an NPSS input file typically follows 
this “recipe” of steps “a” through “f ’. An input file may repeat certain steps in order to analyze many 
cases, the first case being the design point analysis and all other cases representing off-design analysis 
points of interest. 


Thermodynamic Performance Analysis — Design and Off-Design 

Cycle analysis of air-breathing engines is done in two stages. First the engine components are sized 
according to the inputs specified by the user; this is called the design case or design analysis. The design 
case can also be considered a reference condition. Normally only one design case is executed for each 
model run; if multiple cases are being run then while the engine is in design mode each case results in a 
completely different engine. 

The second stage of cycle analysis is off-design analysis. Here an engine whose design variables are 
known has its operating characteristics calculated for a particular flight Mach number, altitude, and 
throttle setting. For this reason off-design analysis is often called performance analysis. Hundreds of off- 
design cases are run to identify the engine operational envelope as well as any unforeseen complications 
resulting from potentially poor design choices. So in design mode the engine physical parameters are set, 
and in off-design mode that engine’s performance is calculated over a range of varying conditions (flight 
speed, altitude, etc.). 


Design Point Analysis 

In order to analyze the performance of an air-breathing engine cycle, a minimum amount of input and 
design conditions must be specified. Table 1 is a listing of the more common engine simulation elements 
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and the minimum necessary input variables for each at the engine design point. Usually turbomachinery 
maps will be used for the performance analysis of compressors and turbines. The two different sets of 
inputs are (not including the map table variables) listed in the bottom half of table 1. Since turbine 
pressure ratio is usually a fallout from design choices rather than a design variable itself, the turbine 
pressure ratio scaler is determined from the desired (unsealed map input) value pressure ratio and the 
actual pressure ratio. See appendix C for more information on turbomachinery map scaling. 


TABLE 1.— CRITICAL INPUT VARIABLES 


Element type 

Variable name(s) 

FlightConditions 

alt, MN, and W 

Inlet 

cRamBase (direct input or subelement) 

Compressor 

see below 

Splitter 

BPR 

FuelStart 

LHV 

Burner 

W fuel/F AR/TtC ombOut 

Turbine 

see below 

Mixer 

stream 1 MN 

Duct 

dPqPbase (direct input or subelement) 

Nozzle 

PsExh or PsExhName, switchType 

Shaft 

Nmech, declared ShaftlnputPorts 

compressor variables (no map) and description 

NpctDes 

percent of physical speed 

eff 

adiabatic efficiency 

PR 

pressure ratio 

| compressor variables (with map) and description 

alpha 

stator/IGV setting or 3-D argument setting 

RlineMap 

the value of R-line to use @ design 

NcDes 

the value of corrected speed to use @ design 

PRdes 

the desired (scaled) value of pressure ratio 

effDes 

the desired (scaled) value of efficiency 

| turbine variables (no map) and description 

eff 

adiabatic efficiency 

PRbase 

pressure ratio 

| turbine variables (with map) and description 

parmGeom 

vane setting or 3-D argument setting 

parmMapDes 

the desired (unsealed) value of PR @ design 

parmNcDes 

the value of corrected speed to use @ design 

parmMap 

the actual (scaled) value of pressure ratio 

effDes 

the desired (scaled) value of efficiency 


At the design point, NPSS automatically balances and matches most of the engine by doing the 
following: a) changing the scale factors on the map corrected flows for all compressors and turbines so 
that each is correctly sized to pass the flow into that component; b) setting the exit area (and if applicable, 
throat area) for all nozzles based on the nozzle entrance flow and back pressure; c) setting the area of 
ducts and flowstations (if Mach number is provided); and d) setting areas for mixer inlet streams by 
matching static pressures and setting mixer outlet total pressure based on conservation of momentum 
(zero net impulse) assuming a constant area mixer. All Elements in NPSS have an option variable called 
“switchDes” which can be set to either “DESIGN” or “OFFDESIGN”; this variable determines whether 
certain variables are input or output (i.e., calculated). 

What normally remains at design is an energy balance for each shaft. The user must add solver 
objects that vary the expansion ratio across each turbine (variable ‘parmMap’) to match the supplied 
turbine work (or torque) to each shaft with the torque required by all compressors and loads also on that 
shaft (variable ‘trqNet’); these solver objects will be added automatically if the autoSolverSetup() 
command is issued at design. Finally, other design conditions can be included if desired. Examples are 
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varying compressor pressure ratio to achieve a certain thrust specific fuel consumption (TSFC) and 
varying engine airflow to meet a certain net thrust. 


Off-Design Point Analysis (Steady State) 

In NPSS, Independent and Dependent objects are used to balance the engine by varying x so that y is 
satisfied, where x is the free variable specified by the Independent and y is the mismatch we want to 
balance specified by the Dependent. Because this is a system of simultaneous equations the number of 
active Independents (variables) must always equal the number of active Dependents (equations) for each 
case, although the number may vary from case to case. 

At off-design there naturally exists a flow mismatch at the entrance to any nozzle or rotating 
component. In addition, there is a shaft power mismatch between the turbine-delivered power and the 
compressor-required power. The engine is balanced at off-design by altering the independent (or free) 
variables of available components to eliminate these discontinuities. 

Table 2 lists the most common free variables available to the user at off-design while table 3 lists 
component errors that need to be resolved for a physically valid solution. So for a two-shaft separate flow 
turbofan each rotating component's free variable is used to balance the flow mismatch at the entrance to 
that particular component and the shaft lpms are used to balance the power delivered/required by each 
shaft. That leaves inlet airflow and bypass ratio to balance the flow mismatch at the two nozzle entrances. 
Any remaining independents are available for the user to control the engine; this is explained further in 
the next section. 


TABLE 2.— AVAILABLE INDEPENDENT VARIABLES AND CORRESPONDING INDEPENDENT OBJECT NAMES 


Element type 

Variable name 

Independent name 

FlightConditions 

W 

ind W 

CompressorMap(subelement) 

RlineMap 

ind RlineMap 

TurbineNeppMap(subelement) 

pannMap 

ind_parmMap 

Shaft 

Nmech 

ind Nmech 

Splitter 

BPR 

ind BPR 

Burner 



(switchBum = WFUEL) 

Wfuel, or 

must be created by user 

(switchBum == TEMPERATURE) 

TtCombOut, or 

must be created by user 

(switchBum == FAR 

FAR 

must be created by user 

FlowStart 

W 

must be created by user 


Pt 

must be created by user 


Tt 

must be created by user 

Mixer 

none — (internal to mixer) 


Nozzle 

none — (fixed area) 



AthCold — (variable area) 

must be created by user 


TABLE 3.— SYSTEM DISCONTINUITY ERRORS AND DEPENDENT NAMES 


Element type 

Error/mismatch 

Dependent name 

CompressorMap 

flow into the compressor 

compressor map flow based on Rline and Nc 

dep_errWc 

TurbineNeppMap 

flow into the turbine 

turbine map flow based on PR and Nc 

dep_errWp 

Nozzle 

flow into the nozzle 

nozzle flow based on area and PR 

dep_Area 

Shaft 

torque/work required for all compressors and loads 
torque/work delivered by all turbines 

integrate Nmech 

Mixer 

( 1 ) stream 1 static pressure 

( 1 ) stream 2 static pressure 

(2) stream 1 impulse + stream 2 impulse 
(2) exit stream impulse * loss factor 

deperrPs 
none — (internal) 
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Off-Design Full Power Operation (Maximum Power) 

In the case of the two-shaft separate flow turbofan mentioned above one free variable remains, burner 
fuel flow, to be used as desired. If burner fuel flow is not explicitly changed or made active in the solver it 
will remain constant. If the engine is operated over a flight envelope this way the turbomachinery (mostly 
the fan and LPC) is forced to odd operating points (e.g., fan at extremely high corrected speed, LPC at 
low surge margin): a way to determine and control fuel flow as engine operating conditions change is 
needed. 

For off-design analysis it is not sufficient to merely balance engine discontinuities, a definition of full 
power operation is required. Examples of full power definitions are maintaining fan corrected speed at 
1 00 percent, maintaining low spool lpm at 1 00 percent, or maintaining a constant high pressure turbine 
entrance temperature. Once the engine full power operating requirements are set, fuel flow can be 
calculated and the analyst can begin to see the effects of design choices on engine performance. 


Off-Design Part Power Operation (Throttle Curves) 

This is a subset of off-design operation. While an engine may operate at maximum power during an 
aircraft takeoff and climb, it is necessary to throttle back the engine during long cruise segments and 
especially during descent. The essence of part power operation is simply a matter of decreasing burner 
exit temperature from the defined maximum at some Mach number and altitude and letting NPSS balance 
the engine. There are several ways of doing this, each with advantages and disadvantages, but all yield the 
same results. Three methods of defining part power operation are explained below. 

The first method is to directly input burner exit conditions. Simply run a series of cases with the 
burner exit temperature/fuel flow/FAR input (which one is input depends on the value of the burner 
variable ‘switchBum’) and decreasing from case to case. This method is simple but it requires prior 
knowledge of the starting (full power) burner temperature/fuel flow/FAR which may or may not be 
constant from case to case. 

The second method is to vary burner exit conditions to yield a given net thrust. To do this an 
Independent and Dependent must be created and made active in the solver. The Independent variable 
comes from the burner Element and depends on the value of the ‘switchBum’ variable. The Dependent 
will equate the engine thrust performance variable to a target thrust value. The target thrust value is 
frequently defined in terms of a percentage of the engine maximum thrust at each flight condition for 
convenience. 

The third method is to vary burner exit conditions to yield a given fan corrected speed (or high 
pressure compressor corrected speed for turbojets). This is similar to the second method; only the 
Dependent conditions are different. 


Common Error Messages in NPSS: What They Mean and How to Fix Them 

NPSS has a large array of warning and error messages to assist the user in eliminating improper 
conditions in a model simulation. When ‘debugging’ a model it is critical that the novice user remedy the 
first error message only! It is very common for a single input error to cascade into a screen-scrolling 
nightmare. The following are some of the more common error messages that occur when constructing a 
model simulation and what causes them. 

ERROR(7 100 1001) in file filename - line n - Unresolved variable ‘xxx’ 

All variables in NPSS must be explicitly declared (as real, int, or string); this error means NPSS has 
encountered a variable which has not been declared. The variable in question could be user-created and 
this error has four likely causes: 
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1) the variable name is misspelled or incorrectly capitalized (this is the likely cause when the 
variable is not user-created) 

2) the variable is simply not declared at all (for user-created variables) 

3) the variable is declared but after NPSS first encounters it being assigned a value 

4) the variable is declared but referenced outside its scope (make sure the scoping name is correctly 
spelled) 

WARNING(71021001) in file filename - line n - in Solver ‘solver’: VERIFY FAILED : Jacobian indep 
and dep lists are not the same length ( #indep = nl, #deps = n2 ) 

This warning indicates the number of active Independents, nl, is not equal to the number of active 
Dependents, n2, at the time of the run() command. An error message will follow this warning. 

ERROR(8 1023999) in file filename - line n - in SecantSolver ‘xxt.smSolver’ : Slope between two points 
went off to infinity: 

For this error message ‘jccc’ is likely to be a rotating component name and the possible cause is that the 
component operating point has moved beyond its map boundaries. Alternatively, if the mechanical speed 
of a shaft (shaft name. N mcch ) connected to a compressor is inadvertently set to zero at design, this error 
will result. 

ERROR(8 103 1001) in file filename - line n - in Assembly function xxxversion was not found in lib 
xxx. dll 

Due to the case sensitivity of NPSS contrasted with the case insensitivity of some computer file systems, 
this error is probably the result of not capitalizing the first letter of an Element type like “compressor” 
instead of “Compressor”. 

ERROR(9 1003 001) in file filename - line n - in MemberFunction 'run': Exception caught in 
SolverExecutive::run : Exception caught in SolverExecutive::convergeContinuous 
Model : Exception thrown from matrix generation step : Exception caught in Jacobian: :generateMatrix 
Exception caught during matrix inversion Exception caught in 

Jacobian: :determineColumnNormalizerSingular Matrix - zero column: Perturbing matrix-independent 
number #nl, independent name, did not effect any of the error terms 

The last part of this long error message is what matters. Each active Independent variable must have an 
effect on at least one active Dependent term. Note that if the dependency is extremely small (but not 
exactly zero) an unconverged solution error will result instead. Assuming that the named Independent can 
significantly affect at least one Dependent then there are two likely causes: 

1) the independent variable is repeatedly being assigned values somewhere other than just the 
solver, as in a preexecute() function 

2) the independent variable value is beyond the range of a table or function whose extrapolation is 
set to “NONE”, resulting in the same tabular value being returned regardless of the independent 
value 

ERROR(9 1003 001) in file filename - line n - in MemberFunction 'run': Exception caught in 
SolverExecutive::run : Exception caught in SolverExecutive::convergeContinuous 
Model : Exception thrown from matrix generation step : Exception caught in Jacobian: :generateMatrix 
Exception caught during matrix inversion Exception caught in Jacobian::determineRowNormalizer 
Singular Matrix - zero row: The error for matrix-row number nl, dependentname, is not affected by any 
of the independent perturbations 

Similar to the previous error, this error occurs when a dependent term remains unchanged by any 
Independent variable change. Two likely causes: 
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1) the Dependent variable terms (“eq_rhs”, “eq_lhs”, or both) are not suitably li nk ed to the active 
Independents 

2) the Dependent variable terms (both “eq_rhs” and “eq_lhs”) are ultimately being assigned to 
constants (e.g., vary x so that 5 equals 7 is impossible) 

ERROR) 9 1003 001) in file filename - line n - in MemberFunction 'run': Exception caught in 
SolverExecutive::run : Exception caught in SolverExecutive::convergeContinuous 
Model : Exception thrown from matrix generation step : Exception caught in 
Jacobian::generateMatrix Exception caught during matrix inversion Exception caught in 
Jacobian: idetennineRowNormalizer Linearly Dependent Matrix: The coefficients for matrix-row 
#nl, dependent jiamel, are a scalar multiple of those for matrix-row #n2, dependent _name2 
This error means that two of the active Dependents are really the same. For example, if x and y are 
variables the two equations x+y = 5 and 2x+2y =10 can not each be used as Dependent conditions since 
they are trivially related. Another possibility is when the same variable can be referenced by different 
names; an example would be a burner connected to a turbine: the same temperature variable can be called 
burner name. Fl O.Tt and turbinename. Fl l.Tt. 

ERROR( 9 1003201) in file filename - line n - in MemberFunction ‘ran’: Exception caught in 
SolverExecutive::run : Solver FAILED TO CONVERGE continuous model in m iterations. 

Unfortunately, this error is as common as it is difficult to determine its exact cause. It occurs when the 
solver is unable to find a solution within the given tolerance and within the number of allowed iterations. 
Making matters worse is this is not a fatal error, allowing model execution to continue and data to be 
output to files where there may be no indication that it is invalid. General things to look for when this 
error is encountered: 

1) typically increasing the allowable number of iterations beyond 100 will not solve the problem 

2) when the design point fails to converge, check the initial guesses and values for the Independent 
variables to make sure they are at least the same order of magnitude as the final values — if 
airflow starts at 0.1 lb m /s, it will take too many perturbations to reach 100 lb m /s 

3) also at design, make sure each Dependent is significantly affected by at least one Independent - 
varying HPC design adiabatic efficiency to achieve a certain HPT design polytropic efficiency is 
an example of a minuscule dependence (HPC efficiency affects HPT pressure ratio which affects 
HPT polytropic efficiency) 

4) for off-design, keep changes in flight Mach and altitude relatively small from case to case — going 
from sea level, static conditions to Mach 0.85, 40000 ft altitude is fine for a simple turbojet but 
far too big a change for complex engine simulations (usually AMN-0. 1 0 and Ah~2000 ft are 
acceptable) 

5) verify that Independent variables are not going beyond table ranges and maps which can result in 
odd extrapolation (or no extrapolation) 

6) occasionally a case will not converge because the model inputs are just right to cause the solver to 
consistently overshoot the solution - if the engine simulation works both at 9990 ft and 10010 ft 
altitudes but not at 1 0000 ft then this is the likely problem 

ERROR(9 104 1001) in file filename - line n - in Assembly Cannot create ‘my object name’ of type 
‘xxx’ 

Similar to the unresolved variable error, this error arises when the user tries to instantiate an object that 
NPSS does not recognize. Likely causes are: 

1) the object type, ‘xxx’, is misspelled like “Nozzel” or “Indepenent” 

2) the object type should be capitalized like “Independent”, not “independent” 
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ERROR(91041001) in file filename - line n - Could not create FlowStation because you haven’t 
specified a thermo package. 

NPSS needs a thermodynamic package in order to create some of its objects; make sure the 
setThermoPackage(“xxx”) command is one of the first things in the input file. 
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Appendix A 

Annotation of Turbojet Example Model 

In the turbojet example, shown schematically in figure 1, the model is composed of eleven Elements, 
seven of which represent the physical components inlet, compressor, burner, turbine, duct, nozzle, and a 
shaft connecting the turbomachinery. The four remaining Elements do not represent physical machinery 
but nevertheless provide important data, set certain conditions, or perform important functions. These four 
Elements are the following: 

FlightConditions — this Element calculates atmospheric properties and starts the fluid flow stream 

FuelStart — this Element specifies fuel properties 

FlowEnd — this Element is required to end the fluid flow stream 

EngPerf — this is an optional Element that calculates engine performance variables 

In general, the minimally required or most important variables for each Element are assigned values. In 
instantiation order, the eleven Elements in the turbojet model are: 

1) FlightConditions Element named “Ambient” 

This Element calculates atmospheric properties based on the two input variables ‘alt’, the engine 
altitude in feet, and ‘MN’, the flight Mach number of the engine. Here the example turbojet is operating at 
sea level (altitude = 0) static (velocity or MN = 0) conditions. The FlightConditions Element also 
specifies the amount of air mass flow into the engine (W =100 lb m /s). 

2) Inlet Element named “Inlet” 

This Element adjusts the fluid pressure according to the inlet stagnation pressure recovery variable, 
‘eRamBase’. In the turbojet example the inlet pressure recovery is 0.980, or a 2 percent total pressure loss 
across the inlet component. 

3) Compressor Element named “HPC” 

This Element is used to represent high pressure compressors, low pressure compressors, and fan 
components in engine models. Because of the large amount of text and information in all but the simplest 
of models, NPSS input is often contained in multiple files for easier study. Any input file can be included 
by the statement #include<my_file_name>, where ‘my file name’ is the name of the file (double 

r-FS 1 r-FS 2 r ~ FS 3 r ~ FS 4 r~FS 5 i ~ FS 7 r ~ FS 9 


i i i i i i i 



Figure 1. — Example turbojet NPSS thermodynamic cycle model. 
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quotes in place of angle brackets may also be used). In the example model, the compressor has a 
CompressorMap Subelement called ‘Smap’ which contains three data tables; this Subelement and 
associated tables are contained in the “turbojetHPC.map” file which is included on line 35 of the model 
file. See appendix C for a discussion of how turbomachinery maps are used in NPSS. In the turbojet 
model the compressor has a pressure ratio of 20.0 and an adiabatic efficiency of 85.0 percent. 

4) FuelStart Element named “Fuelln” 

The purpose of the FuelStart Element is to specify the properties of the fuel that will be added and 
burned in the combustor. The default fuel is JP-4, and the most common input variable is the fuel lower 
heating value, or LHV. The fuel in the example model has an LHV of 1 8400 BTU/lb m . This Element can 
also be considered as starting a flow stream of fuel, although the flow rate is usually set in the Burner 
Element. 

5) Burner Element named “Burner” 

The Burner Element adds fuel to an airflow and bums it, resulting in a temperature rise. The Burner 
calculations depend on the ‘switchBum’ variable: if ‘switchBum’ = “TEMPERATURE”, the Burner exit 
temperature is set equal to the value of the variable ‘TtCombOuf and the necessary fuel-to-air ratio 
(FAR) and fuel flow are calculated. If ‘switchBum’ = “FUEL”, the Burner adds fuel equal to the value of 
the variable ‘Wfuel’ and the corresponding fuel-to-air ratio (FAR) and exit temperature are calculated. If 
‘switchBum’ = “FAR”, the Burner exit fuel-to-air ratio is set equal to the value of the variable ‘FAR’ and 
the corresponding fuel flow (Wfuel) and exit temperature are calculated. This last method is used in the 
turbojet example to set the Burner exit fuel-to-air-ratio to a value of 0.030. Recalling that the input 
airflow is 100 lb m /s, the fuel flow for this model will be 3.0 lb m /s. 

The Burner also has a stagnation pressure drop specified by the variable ‘dPqPBase’ (i.e., AP/P). The 
turbojet example has a stagnation pressure loss of 5 percent across the Burner relative to the Burner 
entrance stagnation pressure. 

6) Turbine Element named “HPT” 

This Element is used to represent both high and low pressure turbine components. Like the 
compressor, the turbine in the example model has Subelements and tabular data contained in an included 
file, “turbojetHPT.map”. See appendix C for more information on turbine maps. The example model 
turbine has an adiabatic efficiency of 92.0 percent. 

7) Duct Element named “Ductl” 

This Element is primarily used to simulate pressure drops from frictional effects in long passages or 
losses from turning the fluid flow. The pressure drop can be set from the input value of ‘dPqPbase’, a 
function of duct entrance Mach number, or a user-specified loss function. The example model has a 
stagnation pressure loss of 1.0 percent. 

8) Nozzle Element named “Nozzle” 

The Nozzle Element represents both convergent and convergent-divergent nozzles as determined by 
the ‘switchType’ variable. Thrust calculations are based on either velocity coefficient or gross thrust 
coefficient as determined by the value of the ‘switchCoef variable. To do the nozzle calculations an 
exhaust or back pressure is required; this can be done by either input of the back pressure itself (variable 
‘PsExh’) or a link to the proper pressure (variable ‘PsExhName’). Normally the nozzle exhausts to the 
atmosphere and this can be specified by setting ‘PsExhName’ equal to the model’s FlightConditions 
Element static pressure variable, as is done in the example model. The example nozzle is convergent 
(‘switchType’ = “CONIC”) with losses based on velocity coefficient (‘switchCoef = “CV” and ‘Cv’ = 
0.98). 
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9) FlowEnd Element named “NozzEnd” 

This Element is used to terminate a flow stream; all flow streams in NPSS must have a beginning 
(either FlowStart or FlightConditions Element) and an end. 

10) Shaft Element named “HPShaft” 

The Shaft Element is used when a mechanical link exists between two or more rotating components. 
Unlike many Elements which are connected by their FluidPorts, this Element connects to others via 
ShaftPorts. By default the Shaft Element has no ShaftPorts; instead its ShaftPorts are declared by the user 
so that it can have any number appropriate to the engine being simulated. The command to instantiate a 
ShaftPort is ShaftlnputPort name, where ‘name’ is a unique name for that port. In the turbojet example 
the shaft element has two ShaftlnputPorts called ‘ShHPC’ and ‘ShHPT’; multiple ports can be 
instantiated in one command using commas to separate the names as in the example. Once the ports have 
been declared, they can be linked; these ports will be connected to the ShaftOutputPorts on the model’s 
compressor and turbine later in the input file. 

The other thing of note in the Shaft Element is the ‘Nmech’ variable which is the shaft rotational 
speed in lpm; the example model has a shaft speed of 1 0000 rpm. 

11) EngPerf Element named “PERF” 

By design, NPSS models can be very different: a model may represent a turboprop, turbofan, just the 
core of an engine, etc. Because of this, a priori equations and variables related to overall engine 
performance for every conceivable model configuration are impossible. However, the EngPerf Element 
calculates some common overall engine performance variables such as net thrust and thrust specific fuel 
consumption. This Element assumes that the model is a thrust-producing configuration (not a turboshaft 
or turboprop) with at least one inlet, at least one burner, and at least one nozzle. In addition, these 
Elements must be the ‘standard’ NPSS release Elements, not user-created Elements. The turbojet example 
can use the EngPerf Element. In the event that the EngPerf Element is not applicable, the user has three 
choices: 

1) omit it and not have any performance variables, 

2) omit it but add their own performance calculations somewhere else, or 

3) modify it to correctly calculate their performance variables 

After all the Elements in the model have been declared they must be linked together forming a 
coherent simulation of a gas turbine engine. The types of ports used in NPSS are FluidPorts, ShaftPorts, 
and ThermalPorts. Ports are linked using the linkPorts command: 
linkPorts( “elementl.outlet port”, “element2.inlet_port”, “linkname” ); 

where ‘link name’ is a user-specified, unique name for the link. By convention, most FluidOutletPorts 
are called “F1_0”, FluidlnputPorts are called “FI I”, and ShaftOutputPorts are “Sh_0”. Finally, 
FluidOutputPorts can only be linked to FluidlnputPorts; likewise, ShaftOutputPorts can only be linked to 
ShaftlnputPorts. In the example model there is a FuelOutputPort on the FuelStart Element connected to 
the FuellnputPort on the Burner Element, and the two declared ShaftlnputPorts on the Shaft Element 
(called “Sh HPC” and “Sh HPT”) are connected to the ShaftOutputPorts on the Compressor and Turbine 
Elements respectively. 

Once all necessary Elements are instantiated, values given to important variables, and all of the 
Element ports are linked, the model simulation can be run. This consists of any necessary changes to the 
solver followed by the run() command. The autoSolverSetup() command makes all Independent and 
Dependent objects inactive then activates any Independent or Dependent whose “autoSetup” variable is 
set to TRUE. 

In the turbojet example, the command setOption(“switchDes”, “DESIGN”) sets the value of every 
variable called “switchDes” in the model to “DESIGN”. A consequence of this is that certain 
Independents and Dependents have their “autoSetup” variable value set to either TRUE or FALSE. In this 
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case the HPT map subelement will have its Independent, “ind_parmMap”, added to the solver and the 
HPShaft Element will have its Dependent, “integrate Nmech”, added to the solver when the 
autoSolverSetup() command is issued because their “autoSetup” value is TRUE. In this example no 
other Independent or Dependent is active during the design case run. 

After the run() command, the cout command is used to display to the screen simple text messages 
and variables related to engine performance. A DataViewer declaration follows, which defines a specific 
type of output format and variables to be shown. In this particular DataViewer, the viewer subtype is a 
PageViewer and the actual name of the viewer is “page”. The name of the file this DataViewer writes to 
is called “turbojet. output”. Declaration of this DataViewer does not send any data to the file; that is done 
when the DataViewer’s display() command is given, as seen on the next line in the example model. 
Notice that the display() command is a function of each individual DataViewer, so it must be given 
proper reference, or scope, by first giving the name of the particular DataViewer whose output is desired. 
In this case the DataViewer name is “page”, hence page.display(). 

Next, two off-design points are run. The command setOption(“switchDes”, “OFFDESIGN”) is 
issued to set each element in off-design mode and update the values of the “autoSetup” variables. The 
autoSolverSetup() command now activates the following solver objects: the HPC map Independent, the 
HPT map Independent, the shaft Independent, the engine airflow Independent, the HPC map Dependent, 
the HPT map Dependent, the shaft Dependent, and the nozzle Dependent. To summarize, the active 
solver objects are: 


active Independents 
“HPC . Smap . indRlineMap” 
“HPT. Smap. ind_parmMap” 
“HPShaft.indNmech” 
“Ambient.ind W” 


active Dependents 
“HPC . Smap . deperrWc” 
“HPT.Smap.deperrWp” 
“HPShaft. integrateNmech” 
“N ozzle. depArea” 


The case is run and data sent to the PageViewer output file. Since the engine operating conditions were 
not changed, the output for this off-design point, case “a”, will exactly match that for the engine’s design 
point. However, for the second off-design point, case “b”, the engine operating conditions are desired at a 
flight Mach number of 1 .00 and a flight altitude of 20000 ft. In addition, the burner exit fiiel-to-air ratio is 
changed to 0.03 14. If the burner FAR is not explicitly changed by the user (or associated with an active 
Independent), it will remain constant. As before, the case is run and the data output. 

As a final example, another off-design point is run with a set of user-defined solver objects active in 
addition to the others. In the previous case, the burner fuel-to-air ratio was set at 0.03 14, yielding a 
combustor exhaust temperature (T 4 ) of approximately 3160 Rankine. Suppose a T 4 equal to 3200 Rankine 
is desired: this could be achieved by trial-and-error input of burner FAR, but creation and use of an 
Independent and Dependent is far superior. Case “c” illustrates this process and runs the case; it can be 
seen from the output that burner FAR has been altered to produce a T 4 of exactly 3200, within solver 
tolerance. Once created, these Independent and Dependents do not have to be specified again. It is 
advisable to create all Independents and Dependents in the same place and use them as necessary for each 
case remembering that they always remain in the same state (whether active or inactive) from case to case 
unless explicitly changed. 
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Appendix B 

Example Turbojet File Listing 


The example turbojet model consists of the following five files: 
turbojet.mdl this file configures the engine and runs the simulation 

turbojetHPC.map this file contains the compressor performance map 

turbojet_HPT.map this file contains the turbine performance map 

turbojet. view_page this file determines the format of the output 

turbojet. output this is the output file generated by running the simulation 

Again, file names can be completely arbitrary but these were chosen based on convention to reflect their 
contents. The “turbojet. view_page” file contains a complex DataViewer specification which uses NPSS 
functions to generate an output format compatible with most model simulations that use the standard 
NPSS Elements. This DataViewer will not be explained in detail, but several lines of interest are 
explained. First, an output stream name must be created with an associated filename where the data will 
be written: 

OutFileStream pageStream { filename = “turbojet. output”; } 

Here the output stream is called ‘pageStream’ and all data in this stream will go to the file 
“turbojet.outpuf ’. Note this is done outside of the actual DataViewer specification, which appears next: 

DataViewer PageViewer page { 

// many lines which list desired output variables in specific groups, and their format 
outStreamHandle = “pageStream”; 

} 

This DataViewer is a PageViewer subtype (VarDumpViewer, CaseColumnViewer, and CaseRowViewer 
are the others) and ‘page’ is the user-given name for this particular viewer. The ‘outStreamHandle’ string 
variable is set equal to the OutFileStream object name previously defined (‘pageStream’), thereby 
associating this viewer’s output with that output stream and the file, “turbojet.output”. Each time this 
DataViewer’s display() function is called, model simulation data is written. In the example model, the 
“turbojet.output” file contains the engine data at four conditions: the DESIGN case and the three OFF- 
DESIGN cases. 


NASA/TM— 2007-2 14690 


15 



turbojet. mdl .txt 


//- 


// 


// 

File Name: 

// 

Date(s) : 

// 

Author : 

// 


// 

Descri pti on : 

// 

// 


//- 



turbojet. mdl 
December 1, 2005 
Scott M. Hones 

simple NPSS turbojet model 


for pedagogical 


use 


setThermoPackage("Hanaf ") ; 

/* 

this is an example of 
multiple line comments 
using the slash and asterisk 

*/ 

Element Fl i ghtCondi ti ons Ambient { 
MN = 0.00; 
alt = 0. ; 

W = 100.0; 

} 


Element Inlet Inlet { 

eRamBase = 0.980; // comment: this is ram recovery 

} 


Element Compressor HPC { 

// use a map for the compressor from the file, "turbojet_HPC.map" 
#include <turbojet_HPC.map> 

S_map . Rl i neMap = 2.0000; 

S_map.NcDes = 1.000; 

S_map.PRdes = 20.00; 

S_map.effDes = 0.8500; 

} 


Element Fuelstart Fuelln { 

LHV = 18400.; // this is in BTU per lb mass 

} 


Element Burner Burner { 
dPqPBase = 0.050; 
switchBurn = "FAR"; 
FAR = 0.0300; 

} 


Element Turbine HPT { 

// use a map for the turbine from the file, "turbojet_HPT.map" 
#include <turbojet_HPT. map> 

S_map.effDes = 0.9200; 

S_map . parmMap = 3.0; 

S_map . parmGeomDes = 1.0; 

S_map . parmMapDes = 3.000; 

S_map . parmNcDes = 100.0; 
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turbojet. mdl .txt 


Element Duct Ductl { 

dPqPbase = 0.010; 

} 


Element Nozzle Nozzle { 
switchType = "CONIC"; 
switchCoef = "CV"; 

Cv = 0.9800; 

PsExhName = "Ambi ent . Ps" ; 

} 


Element FlowEnd NozzEnd { } 


Element Shaft HP_Shaft { 

ShaftlnputPort Sh_HPC, Sh_HPT; 
Nmech = 10000; 

} 


Element EngPerf PERF { 

} 


// 

// 

// 

linkPortsC 
linkPortsC 
1 i nkPorts( 
linkPortsC 
linkPortsC 
li nkPorts( 
linkPortsC 


COMPONENT LINKAGES 


Ambi ent . Fl_0 
Inlet. Fl_0" 
HPC. Fl_0" 
Burner . Fl_0" 
HPT. Fl_0" 
DUCtl. Fl_0" 
Nozzl e . Fl_0" 


Inlet. Fl_l" 
HPC. Fl_l" 
Burner . Fl_l" 
HPT. Fl_l" 
DUCtl. Fl_l" 
Nozzle. Fl_l" 
NozzEnd . Fl_l 


" FS_1" ) 

" FS_2" ) 

" FS_3" ) 

" FS_4" ) 

" FS_5" ) 

" FS_7" ) 

" FS_9" ) 


linkPortsC "Fuel In . Fu_0" 
linkPortsC "HPC.Sh_0" 
linkPortsC "HPT.Sh_0" 


Burner . Fu_I 
HP_Shaft.Sh_HPC 
HP_Shaft . Sh_HPT 


Fuel_l") ; 
HPCworkin" ) ; 
HPTworkout" ) ; 


// RUN THE MODEL DESIGN CASE 


setOption( "switchDes", "DESIGN" ); 
autoSol verSetupO ; 
run() ; 

// WRITE SOME TEXT AND DATA TO THE SCREEN 

cout « "Turbojet Thrust = " « PERF.Fn « " lbs " « endl ; 

cout « "Turbojet Fuel Flow = " « PERF.Wfuel « " lbs/hr " « endl; 

cout « "Turbojet TSFC = " « PERF.TSFC « " lb/hr/lb " « endl; 


// 

// PRINT OUT THE DATA 

// 

#include <turbojet.view_page> 
page . di spl ay() ; 


Page 2 
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turbojet. mdl .txt 


// 

// RUN AN OFF-DESIGN CASE 

// a) AT THE SAME CONDITIONS AS THE DESIGN CASE 


setOption( "switchDes", "OFFDESIGN" ); 
autoSol verSetupO ; 
run() ; 

page . di spl ay() ; 


// 

// RUN AN OFF-DESIGN CASE 

// b) AT A DIFFERENT ALTITUDE AND MACH NUMBER 


Ambient. alt = 20000.; 
Ambient. MN = 1.00; 
Burner. FAR = 0.0314; 
run() ; 

page . di spl ay() ; 


// 

// 

// 

// 


RUN AN OFF-DESIGN CASE 

c) USING NEWLY-CREATED SOLVER OBJECTS AND MAKING THEM ACTIVE 


Independent Burner_FAR { 
varName = "Burner . FAR" ; 

} 


// <-- the Independent's name is "Burner_FAR" 
// <-- the variable that this Indep. can vary 
// this Independent has now been CREATED 


Dependent Target_T4 { // <-- the Dependent's name is "Target_T4" 

eq_rhs = "Burner . Fl_O.Tt" ; // <-- right-hand-side of the Dep. condition 

eq_lhs = "3200.0"; // <-- left-hand-side of the Dep. condition 

} // this Dependent has now been CREATED 


sol ver . addlndependentC "Burner_FAR" ); // this Independent is now ACTIVE 

sol ver . addDependent( "Target_T4" ); // this Dependent is now ACTIVE 

run() ; 

page . di spl ay() ; 


Page 3 
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tu r bo j et_HPC . map . txt 


// File Name: turbojet_HPC.map 

// This data Map is generated from an existing NEPP map. 

// The original NEPP map tables were: 

// 3001 BASELINE HP COMPRESSOR FLOW VS. R, SPEED, AND ANGL 

// 3002 BASELINE HP COMPRESSOR EFF VS. R, SPEED, AND ANGL 

// 3003 BASELINE HP COMPRESSOR PR VS . R, SPEED, AND ANGL 

// 

// 


// 

// | Declaration of a new subelement instance compressorMap 

// j that is of the type: CompressorMap 

// 


Subelement CompressorMap S_map { 
al pha = 0.0; 
a_alpha = 0.0; 
s_alpha = 1.0; 

RlineMap = 2.000; 

NcDes = 1.000; 

PRdes = 23.00; 
effDes = 0.8520; 

ReDes = 5000000.0; 

RtMap = 0.0685512; 
gamtMap = 1.40052; 

SMWmap = 15.0; 

SMNmap = 18.0; 


// 

// 

// 

// 


Declaration of a new subelement instance compressorEffi ci encyMap 
that is of the type: CompressorEffi ci encyMap 


Subelement CompressorEffi ci encyMap S_eff { 
RlineStall= 1.0; 


Table TB_Wc(real alphaMap, real NcorrMap, real RlineMap) { 
alphaMap= 0.000 { 


NcorrMap= 0.500 { 


RlineMap = { 

1.000, 

1.200, 

1.400, 

1.600, 

1.800, 

2.000, 

2.200, 

2.400, 

2.600, 

2.800, 

3.000 } 


WcorrMap = { 

22.7411, 

24.0487, 

25.1548, 

26.0615, 

26.7738, 

27.2992, 

} 

NcorrMap= 0.600 

27.6470, 

27.8286, 

27.8634, 

27.8634, 

27.8634 } 

{ 






RlineMap = *; 
WcorrMap = { 

31.7548, 

33.1181, 

34.2670, 

35.2054, 

35.9397, 

36.4783, 

} 

NcorrMap= 0.700 

36.8308, 

37.0085, 

37.0362, 

37.0362, 

37.0362 } 


{ 






RlineMap = *; 
WcorrMap = { 

46.1066, 

47.4088, 

48.5066, 

49.4046, 

50.1096, 

50.6291, 

} 

NcorrMap= 0.750 

50.9717, 

51.1469, 

51.1757, 

51.1757, 

51.1757 } 


{ 






RlineMap = *; 







WcorrMap = { 

56.7268, 

58.0480, 

59.1608, 

60.0704, 

60.7837, 

61.3084, 

61.6527, 

61.8260, 

61.8517, 

61.8517, 

61.8517 } 



Page 1 
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tu r bo j et_HPC . map . txt 


NcorrMap= 0.800 { 
RlineMap = *; 

WcorrMap = { 70.1448, 

75.2269, 

} 

NcorrMap= 0.850 { 
RlineMap = *; 

WcorrMap = { 89.3764, 

95.1897, 

} 

NcorrMap= 0.900 { 
RlineMap = *; 

WcorrMap = { 118.0620, 
125.3385, 

} 

NcorrMap= 0.925 { 
RlineMap = *; 

WcorrMap = { 138.5093, 
146.8174, 

} 

NcorrMap= 0.950 { 
RlineMap = *; 

WcorrMap = { 160.6243, 
167.5334, 

} 

NcorrMap= 0.975 { 
RlineMap = *; 

WcorrMap = { 181.7993, 
187.9029, 

} 

NcorrMap= 1.000 { 
RlineMap = *; 

WcorrMap = { 202.6315, 
206.2702, 

} 

NcorrMap= 1.025 { 
RlineMap = *; 

WcorrMap = { 209.9986, 
212.3516, 

} 

NcorrMap= 1.050 { 
RlineMap = *; 

WcorrMap = { 216.6847, 
218.1106, 

} 


71.5163, 

72.6688, 

73.6088, 

74.3429, 

74.8795, 

75.3943, 

75.4134, 

75.4134, 

75.4134 

} 

90.9746, 

92.3098, 

93.3900, 

94.2232, 

94.8199, 

95.3442, 

95.3504, 

95.3504, 

95.3504 

} 

120.1207, 

121.8253, 

123.1867, 

124.2166, 

124.9292, 


.4609, 125.4609, 125.4609, 125.4609 } 


,8966, 142.8639, 
,9192, 146.9192, 


,5676, 164.1805, 
,6563, 167.6563, 


.4993, 184.9150, 
.0273, 188.0271, 


.5858, 204.3958, 
.4145, 206.4418, 


.5917, 211.1029, 
.4735, 212.5220, 


.0279, 217.3287, 
.2041, 218.2586, 


144.4238, 145.5916, 146.3836, 
146.9192, 146.9192 } 


165.4722, 166.4536, 167.1370, 
167.6563, 167.6563 } 


186.0545, 186.9260, 187.5389, 
188.0271, 188.0271 } 


205.0661, 205.5998, 206.0000, 
206.4418, 206.4418 } 


211.5321, 211.8825, 212.1554, 
212.5227, 212.5227 } 


217.5860, 217.8015, 217.9767, 
218.2739, 218.2739 } 


alphaMap= 90.000 { 
NcorrMap= 0.500 { 

RlineMap = { 1.000, 1 
2.400, 2.600, 2 
WcorrMap = { 45.5139 

53.4791 

} 

NcorrMap= 0.600 { 
RlineMap = *; 

WcorrMap = { 61.9648 

70.0509 

} 

NcorrMap= 0.700 { 
RlineMap = *; 

WcorrMap = { 84.8030 

92.2351 


. 200 , 

.800, 

, 47. 

, 53. 


1.400, 1.600, 1. 
3.000 } 

7923, 49.6744, 

5986, 53.5986, 


,2862, 66.1998, 

,1629, 70.1629, 


,9149, 88.6608, 

,3561, 92.3561, 
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800, 2.000, 2.200, 

51.1682, 52.2877, 53.0508, 

53.5986, 53.5986 } 


67.7166, 68.8513, 69.6225, 

70.1629, 70.1629 } 


90.0516, 91.1005, 91.8227, 

92.3561, 92.3561 } 
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} 


NcorrMap= 0. 
RlineMap = 
WcorrMap = 


750 { 

: {’l00. 
107. 


NcorrMap= 0. 
RlineMap = 
WcorrMap = 


800 { 

: { ' 116 . 
124. 


NcorrMap= 0 
RlineMap = 
WcorrMap = 


850 { 

= { ' 135 . 
143. 


NcorrMap= 0 
RlineMap = 
WcorrMap = 


900 { 

= { ' 157 . 
165. 


NcorrMap= 0. 
RlineMap = 
WcorrMap = 


925 { 

: { ' 167 . 
177. 


NcorrMap= 0. 
RlineMap = 
WcorrMap = 


950 { 

: { ' 180 . 
187. 


NcorrMap= 0 
RlineMap = 
WcorrMap = 


975 { 

= { ' 191 . 
197. 


NcorrMap= 1 
RlineMap = 
WcorrMap = 


000 { 

= {’ 202 . 
206. 


NcorrMap= 1. 
RlineMap = 
WcorrMap = 


025 { 

: {’209. 
212 . 


NcorrMap= 1. 
RlineMap = 
WcorrMap = 


050 { 

: {’216. 
218. 


102.1580, 

107.5439, 


118.8724, 

124.2228, 


137.9187, 

143.6591, 


159.6273, 

166.0495, 


170.4948, 

177.2440, 


182.6532, 

188.0393, 


192.9175, 

197.5237, 


203.5858, 

206.4145, 


210.5917, 

212.4735, 


217.0279, 

218.2041, 


103.8778, 105.2502, 
107.5439, 107.5439, 


120.5804, 121.9443, 
124.2228, 124.2228, 


139.7752, 141.2518, 
143.6591, 143.6591, 


161.7610, 163.4443, 
166.0495, 166.0495, 


172.7683, 174.5537, 
177.2440, 177.2440, 


184.4021, 185.7931, 
188.0393, 188.0393, 


194.3772, 195.5477, 
197.5237, 197.5237, 


204.3958, 205.0661, 
206.4418, 206.4418, 


211.1029, 211.5321, 
212.5220, 212.5227, 


217.3287, 217.5860, 
218.2586, 218.2739, 


106.2882, 107.0060, 
107.5439 } 


122.9763, 123.6901, 
124.2228 } 


142.3627, 143.1228, 
143.6591 } 


164.6937, 165.5281, 
166.0495 } 


175.8698, 176.7370, 
177.2440 } 


186.8382, 187.5511, 
188.0393 } 


196.4376, 197.0574, 
197.5237 } 


205.5998, 206.0000, 
206.4418 } 


211.8825, 212.1554, 
212.5227 } 


217.8015, 217.9767, 
218.2739 } 


Table TB_eff(real alphaMap, real NcorrMap, real RlineMap) { 
alphaMap= 0.000 { 

NcorrMap= 0.500 { 

RlineMap = { 1.000, 1.200, 1.400, 1.600, 1.800, 2.000, 2.200, 

2.400, 2.600, 2.800, 3.000 } 

effAdiabMap = { 0.6753, 0.6913, 0.7016, 0.7050, 0.7004, 0.6864, 0.6570, 
0.6044, 0.5236, 0.4075, 0.2467 } 

} 

NcorrMap= 0.600 { 

Page 3 


NASA/TM— 2007-2 14690 


21 



tu r bo j et_HPC . map . txt 


RlineMap = *; 
effAdiabMap = 

} 

NcorrMap= 0.700 { 
RlineMap = *; 
effAdiabMap = { 

} 

NcorrMap= 0.750 { 
RlineMap = *; 
effAdiabMap = { 

} 

NcorrMap= 0.800 { 
RlineMap = *; 
effAdiabMap = { 

} 

NcorrMap= 0.850 { 
RlineMap = *; 
effAdiabMap = { 

} 

NcorrMap= 0.900 { 
RlineMap = *; 
effAdiabMap = { 

} 

NcorrMap= 0.925 { 
RlineMap = *; 
effAdiabMap = { 

} 

NcorrMap= 0.950 { 
RlineMap = *; 
effAdiabMap = { 

} 

NcorrMap= 0.975 { 
RlineMap = *; 
effAdiabMap = { 

} 

NcorrMap= 1.000 { 
RlineMap = *; 
effAdiabMap = { 

} 

NcorrMap= 1.025 { 
RlineMap = *; 
effAdiabMap = { 

} 

NcorrMap= 1.050 { 
RlineMap = *; 
effAdiabMap = { 

} 

} 

alphaMap= 90.000 { 

NcorrMap= 0.500 { 
RlineMap = { 


{ 0.6953, 0.7094, 
0.6378, 0.5717, 


0.7248, 0.7359, 
0.6838, 0.6366, 


0.7427, 0.7533, 
0.7108, 0.6703, 


0.7634, 0.7736, 
0.7394, 0.7041, 


0.7891, 0.8008, 
0.7754, 0.7417, 


0.8139, 0.8280, 
0.8117, 0.7779, 


0.8206, 0.8356, 
0.8229, 0.7892, 


0.8403, 0.8512, 
0.8415, 0.8179, 


0.8408, 0.8492, 
0.8394, 0.8209, 


0.8470, 0.8505, 
0.8418, 0.8324, 


0.8350, 0.8364, 
0.8275, 0.8217, 


0.8202, 0.8203, 
0.8124, 0.8088, 


1 . 000 , 1 . 200 , 


0.7184, 0.7214, 0.7176, 0.7058, 0.6812, 
0.4783, 0.3512 } 


0.7429, 0.7452, 0.7424, 0.7335, 0.7154, 
0.5713, 0.4848 } 


0.7600, 0.7627, 0.7606, 0.7533, 0.7379, 
0.6147, 0.5414 } 


0.7804, 0.7834, 0.7822, 0.7762, 0.7630, 
0.6556, 0.5920 } 


0.8090, 0.8134, 0.8136, 0.8092, 0.7974, 
0.6950, 0.6335 } 


0.8385, 0.8449, 0.8469, 0.8439, 0.8333, 
0.7303, 0.6671 } 


0.8469, 0.8541, 0.8567, 0.8544, 0.8442, 
0.7416, 0.6783 } 


0.8593, 0.8643, 0.8660, 0.8641, 0.8566, 
0.7852, 0.7423 } 


0.8552, 0.8588, 0.8597, 0.8578, 0.8516, 
0.7954, 0.7624 } 


0.8529, 0.8539, 0.8536, 0.8520, 0.8483, 
0.8200, 0.8043 } 


0.8371, 0.8370, 0.8362, 0.8346, 0.8318, 
0.8141, 0.8049 } 


0.8201, 0.8195, 0.8185, 0.8171, 0.8152, 
0.8045, 0.7992 } 


1.400, 1.600, 1.800, 2.000, 2.200, 
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2.400, 2.600, 

effAdiabMap = { 0.6993, 0.7223, 
0.6580, 0.5734, 

} 

NcorrMap= 0.600 { 

RlineMap = *; 

effAdiabMap = { 0.7294, 0.7502, 
0.7036, 0.6363, 

} 

NcorrMap= 0.700 { 

RlineMap = *; 

effAdiabMap = { 0.7720, 0.7887, 
0.7592, 0.7121, 

} 

NcorrMap= 0.750 { 

RlineMap = *; 

effAdiabMap = { 0.7928, 0.8080, 
0.7855, 0.7457, 

} 

NcorrMap= 0.800 { 

RlineMap = *; 

effAdiabMap = { 0.8125, 0.8267, 
0.8099, 0.7757, 

} 

NcorrMap= 0.850 { 

RlineMap = *; 

effAdiabMap = { 0.8214, 0.8356, 
0.8224, 0.7901, 

} 

NcorrMap= 0.900 { 

RlineMap = *; 

effAdiabMap = { 0.8276, 0.8430, 
0.8332, 0.8005, 

} 

NcorrMap= 0.925 { 

RlineMap = *; 

effAdiabMap = { 0.8245, 0.8403, 
0.8306, 0.7975, 

} 

NcorrMap= 0.950 { 

RlineMap = *; 

effAdiabMap = { 0.8341, 0.8451, 
0.8352, 0.8117, 

} 

NcorrMap= 0.975 { 

RlineMap = *; 

effAdiabMap = { 0.8379, 0.8463, 
0.8365, 0.8179, 

} 

NcorrMap= 1.000 { 

RlineMap = *; 

effAdiabMap = { 0.8470, 0.8505, 
0.8418, 0.8324, 

} 

NcorrMap= 1.025 { 

RlineMap = *; 

effAdiabMap = { 0.8350, 0.8364, 
0.8275, 0.8217, 

} 

NcorrMap= 1.050 { 

RlineMap = *; 

effAdiabMap = { 0.8202, 0.8203, 
0.8124, 0.8088, 


NASA/TM— 2007-2 14690 


2.800, 3.000 } 

0.7389, 0.7480, 0.7483, 0.7382, 0.7110, 
0.4498, 0.2768 } 


0.7653, 0.7739, 0.7751, 0.7676, 0.7460, 
0.5394, 0.4062 } 


0.8009, 0.8080, 0.8094, 0.8044, 0.7892, 
0.6454, 0.5558 } 


0.8192, 0.8259, 0.8276, 0.8237, 0.8109, 
0.6896, 0.6147 } 


0.8372, 0.8436, 0.8456, 0.8426, 0.8318, 
0.7275, 0.6635 } 


0.8463, 0.8529, 0.8553, 0.8528, 0.8429, 
0.7446, 0.6843 } 


0.8546, 0.8622, 0.8652, 0.8633, 0.8537, 
0.7544, 0.6931 } 


0.8521, 0.8599, 0.8630, 0.8611, 0.8514, 
0.7507, 0.6887 } 


0.8531, 0.8581, 0.8598, 0.8578, 0.8504, 
0.7790, 0.7364 } 


0.8523, 0.8559, 0.8568, 0.8549, 0.8486, 
0.7925, 0.7596 } 


0.8529, 0.8539, 0.8536, 0.8520, 0.8483, 
0.8200, 0.8043 } 


0.8371, 0.8370, 0.8362, 0.8346, 0.8318, 
0.8141, 0.8049 } 


0.8201, 0.8195, 0.8185, 0.8171, 0.8152, 
0.8045, 0.7992 } 
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} 

} 

Table TB_PR(real alphaMap, real NcorrMap, real RlineMap) { 


alphaMap= 0.000 { 
NcorrMap= 0 . 500 
RlineMap = { 

PratioMap = { 

} 

NcorrMap= 0.600 
RlineMap = *; 
PratioMap = { 

} 

NcorrMap= 0.700 
RlineMap = *; 
PratioMap = { 

} 

NcorrMap= 0.750 
RlineMap = *; 
PratioMap = { 

} 

NcorrMap= 0.800 
RlineMap = *; 
PratioMap = { 

} 

NcorrMap= 0.850 
RlineMap = *; 
PratioMap = { 

} 

NcorrMap= 0.900 
RlineMap = *; 
PratioMap = { 

} 

NcorrMap= 0.925 
RlineMap = *; 
PratioMap = { 

} 

NcorrMap= 0.950 
RlineMap = *; 
PratioMap = { 

} 

NcorrMap= 0.975 
RlineMap = *; 
PratioMap = { 

} 

NcorrMap= 1.000 
RlineMap = *; 
PratioMap = { 

} 

NcorrMap= 1.025 


{ 


1 . 000 , 

1 . 200 , 

2 . 200 , 

2 . 400 , 

2 . 4769 , 

2 . 4288 , 

1 . 9284 , 

1 . 7711 , 

{ 

3 . 4633 , 

3 . 3778 , 

2 . 5679 , 

2 . 3253 , 

{ 

5 . 0821 , 

4 . 9375 , 

3 . 7106 , 

3 . 3602 , 

{ 

6 . 3490 , 

6 . 1658 , 

4 . 6377 , 

4 . 2042 , 

{ 

8 . 0021 , 

7 . 7686 , 

5 . 8504 , 

5 . 3097 , 

{ 

10 . 4899 , 

10 . 1976 , 

7 . 6539 , 

6 . 9201 , 

{ 

14 . 4564 , 

14 . 0970 , 

10 . 5377 , 

9 . 4621 , 

{ 

17 . 4426 , 

17 . 0500 , 

12 . 7692 , 

11 . 4347 , 

{ 

20 . 7403 , 

20 . 2486 , 

15 . 7626 , 

14 . 4263 , 

{ 

23 . 8298 , 

23 . 2601 , 

18 . 5212 , 

17 . 1524 , 

{ 

26 . 6962 , 

26 . 0933 , 

22 . 0495 , 

20 . 9930 , 


{ 


1 . 400 , 

2 . 600 , 

2 . 3620 , 

1 . 5958 , 

1 . 600 , 

2 . 800 , 

2 . 2778 , 

1 . 4083 , 

3 . 2643 , 

2 . 0595 , 

3 . 1248 , 

1 . 7802 , 

4 . 7554 , 

2 . 9800 , 

4 . 5391 , 

2 . 5826 , 

5 . 9371 , 

3 . 7342 , 

5 . 6667 , 

3 . 2431 , 

7 . 4792 , 

4 . 7237 , 

7 . 1388 , 

4 . 1114 , 

9 . 8249 , 

6 . 1229 , 

9 . 3786 , 

5 . 2899 , 

13 . 6074 , 

8 . 2878 , 

12 . 9977 , 

7 . 0614 , 

16 . 4870 , 

9 . 9718 , 

15 . 7661 , 

8 . 4425 , 

19 . 6093 , 

12 . 9562 , 

18 . 8329 , 

11 . 3983 , 

22 . 5536 , 

15 . 6466 , 

21 . 7200 , 

14 . 0424 , 

25 . 4175 , 

19 . 8436 , 

24 . 6733 , 

18 . 6163 , 
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1 . 800 , 

2 . 000 , 

3.000 } 


2 . 1774 , 

2 . 0627 , 

1.2146 } 


2 . 9619 , 

2 . 7787 , 

1.4973 } 


4 . 2923 , 

4 . 0194 , 

2.1813 } 


5 . 3594 , 

5 . 0204 , 

2.7467 } 


6 . 7532 , 

6 . 3287 , 

3.4919 } 


8 . 8669 , 

8 . 2989 , 

4.4495 } 



12 . 2808 , 11 . 4715 , 
5.8306 } 


14 . 9034 , 13 . 9183 , 
6.9104 } 


17 . 9324 , 16 . 9227 , 
9.8011 } 


20 . 7705 , 19 . 7178 , 
12.3810 } 


23 . 8656 , 22 . 9999 , 
17.3267 } 
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RlineMap = *; 
PratioMap = { 

} 

NcorrMap= 1.050 
RlineMap = *; 
PratioMap = { 

} 


27.6439, 

23.6033, 

{ 

28.4663, 

25.1640, 


27.0687, 26.4522, 25.7969, 25.1052, 
22.7614, 21.8600, 20.9058, 19.9057 } 


27.9667, 27.4460, 26.9054, 26.3460, 
24.5225, 23.8472, 23.1399, 22.4038 } 


} 

alphaMap= 90.000 { 
NcorrMap= 0.500 { 


RlineMap = { 

1.000, 

1.200, 

1.400, 

1.600, 

1.800, 

2.200, 

2.400, 

2.600, 

2.800, 

3.000 } 

PratioMap = { 

5.5200, 

5.4081, 

5.2155, 

4.9486, 

4.6162, 

} 

NcorrMap= 0.600 

3.7722, 

3.2417, 

2.6710, 

2.0954, 

1.5493 } 

{ 





RlineMap = *; 
PratioMap = { 

7.5640, 

7.3974, 

7.1268, 

6.7607, 

6.3104, 

} 

NcorrMap= 0.700 

5.1790, 

4.4722, 

3.7116, 

2.9418, 

2.2066 } 

{ 





RlineMap = *; 
PratioMap = { 

10.2643, 

10.0087, 

9.6410, 

9.1703, 

8.6086, 

} 

NcorrMap= 0.750 

7.2293, 

6.3765, 

5.4525, 

4.5013, 

3.5672 } 

{ 





RlineMap = *; 
PratioMap = { 

12.0520, 

11.7493, 

11.3275, 

10.7965, 

10.1683, 

} 

NcorrMap= 0.800 

8.6354, 

7.6892, 

6.6601, 

5.5926, 

4.5321 } 

{ 





RlineMap = *; 
PratioMap = { 

14.2753, 

13.9203, 

13.4356, 

12.8317, 

12.1210, 

} 

NcorrMap= 0.850 

10.3925, 

9.3260, 

8.1618, 

6.9465, 

5.7278 } 

{ 





RlineMap = *; 
PratioMap = { 

16.9112, 

16.5159, 

15.9649, 

15.2698, 

14.4454, 

} 

NcorrMap= 0.900 

12.4221, 

11.1638, 

9.7852, 

8.3414, 

6.8896 } 

{ 





RlineMap = *; 
PratioMap = { 

20.7368, 

20.3081, 

19.6692, 

18.8350, 

17.8246, 

} 

NcorrMap= 0.925 

15.2936, 

13.6955, 

11.9366, 

10.0936, 

8.2454 } 

{ 





RlineMap = *; 
PratioMap = { 

22.5680, 

22.1257, 

21.4443, 

20.5401, 

19.4347, 

} 

NcorrMap= 0.950 

16.6421, 

14.8683, 

12.9140, 

10.8676, 

8.8203 } 

{ 





RlineMap = *; 
PratioMap = { 

24.1283, 

23.5795, 

22.8503, 

21.9531, 

20.9035, 

} 

NcorrMap= 0.975 

18.3522, 

16.7697, 

15.0246, 

13.1743, 

11.2785 } 

{ 





RlineMap = *; 






PratioMap = { 

25.7180, 

25.1128, 

24.3564, 

23.4593, 

22.4338, 

19.9953, 

18.5065, 

16.8666, 

15.1186, 

13.3083 } 


} 

NcorrMap= 1.000 { 
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24.3798, 


25.7690, 


2 . 000 , 

4.2292, 


5.7900, 


7.9699, 


9.4573, 


11.3187, 


13.5091, 


16.6613, 


18.1546, 


19.7197, 


21.2942, 
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RlineMap : 
Prati oMap 

} 

NcorrMap= 1 
RlineMap : 
Prati oMap 

} 

NcorrMap= 1 
RlineMap : 
Prati oMap 

} 


= { 26 . 6962 , 26 . 0933 , 25 . 4175 , 24 . 6733 , 23 . 8656 , 22 . 9999 , 
22 . 0495 , 20 . 9930 , 19 . 8436 , 18 . 6163 , 17.3267 } 

025 { 

= { 27 . 6439 , 27 . 0687 , 26 . 4522 , 25 . 7969 , 25 . 1052 , 24 . 3798 , 
23 . 6033 , 22 . 7614 , 21 . 8600 , 20 . 9058 , 19.9057 } 

050 { 

= { 28 . 4663 , 27 . 9667 , 27 . 4460 , 26 . 9054 , 26 . 3460 , 25 . 7690 , 
25 . 1640 , 24 . 5225 , 23 . 8472 , 23 . 1399 , 22.4038 } 
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// 

// File Name: turbojet_HPT.map 

// This data Map is generated from an existing NEPP map. 

// The original NEPP map tables were: 

// 1269 BASELINE HP TURBINE EFFICIENCY VS. PR, RPM, AND AREA 

// 1270 BASELINE HP TURBINE FLOW FUNCTION VS. PR, RPM, AND AREA 

// 

// 


// 

// 

// 

// 


Declaration of a new SubElement instance called S_map 
that is of the type: S_map 


Subelement Turbi neNeppMap S_map { 

s_parmMap = 1.0; 
a_parmMap = 0.0; 
s_parmGeom = 1.0; 
a_parmGeom = 0.0; 
s_rni = 1.0; 
a_RNI =0.0; 
effDes = 0.9328; 
parmGeomDes= 1.0; 
parmMapDes = 5.000; 
parmNcDes = 100.0; 
parmGeom = 1.0; 
parmMap = 5.000; 


// 

// | Declaration of a new SubElement instance called Turbi neEffi ci encyMap 

// j that is of the type: Turbi neEffi ci encyMap 

// 

Subelement Turbi neEffi ci encyMap S_eff { 


Table TB_eff(real PgeomDesign, real NcDes , real PRdes) { 


PgeomDesign = 1.0 { 
NcDes= 60.000 { 


PRdes = { 3.000, 

3.250, 

3.500, 

3.750, 

4.000, 

4.250, 

4.500, 

4.750, 

5.000, 

5.250, 

5.500, 

5.750, 

6.000, 

6.250, 

6.500, 

6.750, 

7.000, 

7.250, 

7.500, 

8.000 } 


effMap = { 0.8460, 

0.8405, 

0.8349, 

0.8296, 

0.8249, 

0.8206, 

0.8165, 

0.8127, 

0.8092, 

0.8060, 

0.8030, 

0.8002, 

0.7976, 

0.7953, 

0.7931, 

} 

NcDes= 70.000 { 

0.7911, 

0.7892, 

0.7875, 

0.7858, 

0.7826 } 


PRdes = *; 
effMap = { 0.8879, 

0.8842, 

0.8804, 

0.8769, 

0.8735, 

0.8701, 

0.8670, 

0.8640, 

0.8614, 

0.8590, 

0.8568, 

0.8548, 

0.8529, 

0.8511, 

0.8495, 

} 

NcDes= 80.000 { 

0.8479, 

0.8460, 

0.8436, 

0.8414, 

0.8373 } 


PRdes = *; 
effMap = { 0.9125, 

0.9111, 

0.9096, 

0.9078, 

0.9055, 

0.9034, 

0.9014, 

0.8995, 

0.8979, 

0.8964, 

0.8950, 

0.8936, 

0.8924, 

0.8903, 

0.8877, 

0.8853, 

0.8830, 

0.8808, 

0.8787, 

0.8749 } 


} 

NcDes= 90.000 { 
PRdes = *; 
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tu r ben et_HPT . map . txt 


effMap = { 0.9228, 

0.9242, 

0.9250, 

0.9217, 

0.9210, 

0.9203, 

0.9113, 

} 

NcDes= 100.000 { 

0.9091, 

0.9070, 

PRrlp ^ — « ■ 

effMap = { 0.9215, 

0.9258, 

0.9289, 

0.9326, 

0.9328, 

0.9329, 

0.9245, 

} 

NcDes= 110.000 { 

0.9225, 

0.9206, 

PRdes = *; 
effMap = { 0.9106, 

0.9176, 

0.9232, 

0.9340, 

0.9351, 

0.9361, 

0.9293, 

} 

0.9276, 

0.9259, 

} 

PgeomDesign = 2.0 { 



NcDes= 60.000 { 



PRdes = { 3.000, 

3.250, 

3.500, 

4.750, 

5.000, 

5.250, 

6.500, 

6.750, 

7.000, 

effMap = { 0.8460, 

0.8405, 

0.8349, 

0.8127, 

0.8092, 

0.8060, 

0.7931, 

} 

NcDes= 70.000 { 

0.7911, 

0.7892, 

PRdes = *; 
effMap = { 0.8879, 

0.8842, 

0.8804, 

0.8640, 

0.8614, 

0.8590, 

0.8495, 

} 

NcDes= 80.000 { 

0.8479, 

0.8460, 

PRdes = *; 
effMap = { 0.9125, 

0.9111, 

0.9096, 

0.8995, 

0.8979, 

0.8964, 

0.8877, 

} 

NcDes= 90.000 { 

0.8853, 

0.8830, 

PRdes = *; 
effMap = { 0.9228, 

0.9242, 

0.9250, 

0.9217, 

0.9210, 

0.9203, 

0.9113, 

} 

NcDes= 100.000 { 

0.9091, 

0.9070, 

PRflp^ — " ■ 

effMap = { 0.9215, 

0.9258, 

0.9289, 

0.9326, 

0.9328, 

0.9329, 

0.9245, 

} 

NcDes= 110.000 { 

0.9225, 

0.9206, 

PRdes = *; 
effMap = { 0.9106, 

0.9176, 

0.9232, 

0.9340, 

0.9351, 

0.9361, 

0.9293, 

} 

} 

PgeomDesi gn . i nterp = " 

0.9276, 

0.9259, 

linear" ; 


PgeomDesi gn . extrap = " 

none" ; 


NcDes.i nterp = "linear 
NcDes. extrap = "none" 

J 

J 

Page 
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0.9247, 

0.9240, 

0.9232, 

0.9224, 

0.9197, 

0.9188, 

0.9162, 

0.9137, 

0.9050, 

0.9031, 

0.8980 } 


0.9304, 

0.9313, 

0.9319, 

0.9323, 

0.9330, 

0.9311, 

0.9288, 

0.9266, 

0.9188, 

0.9161, 

0.9107 } 


0.9267, 

0.9292, 

0.9312, 

0.9327, 

0.9369, 

0.9349, 

0.9329, 

0.9311, 

0.9239, 

0.9212, 

0.9161 } 



3.750, 

4.000, 

4.250, 

4.500 

5.500, 

5.750, 

6.000, 

6.250 

7.250, 

7.500, 

8.000 } 


0.8296, 

0.8249, 

0.8206, 

0.8165 

0.8030, 

0.8002, 

0.7976, 

0.7953 

0.7875, 

0.7858, 

0.7826 } 



0.8769, 

0.8735, 

0.8701, 

0.8670 

0.8568, 

0.8548, 

0.8529, 

0.8511 

0.8436, 

0.8414, 

0.8373 } 


0.9078, 

0.9055, 

0.9034, 

0.9014 

0.8950, 

0.8936, 

0.8924, 

0.8903 

0.8808, 

0.8787, 

0.8749 } 


0.9247, 

0.9240, 

0.9232, 

0.9224 

0.9197, 

0.9188, 

0.9162, 

0.9137 

0.9050, 

0.9031, 

0.8980 } 


0.9304, 

0.9313, 

0.9319, 

0.9323 

0.9330, 

0.9311, 

0.9288, 

0.9266 

0.9188, 

0.9161, 

0.9107 } 


0.9267, 

0.9292, 

0.9312, 

0.9327 

0.9369, 

0.9349, 

0.9329, 

0.9311 

0.9239, 

0.9212, 

0.9161 } 




tu r bo j et_HPT . map . txt 


PRdes.interp = "linear" ; 
PRdes. extrap = "none" ; 


Table TB_Wp(real PgeomDesign, real NcDes , real PRdes) { 

PgeomDesign = 1.0 { 

NcDes= 60.000 { 


PRdes 

{ 3 . 000 , 

3 . 250 , 

3 . 500 , 

3 . 750 , 

4 . 000 , 

4 . 250 , 

4 . 500 , 


4 . 750 , 

5 . 000 , 

5 . 250 , 

5 . 500 , 

5 . 750 , 

6 . 000 , 

6 . 250 , 


6 . 500 , 

6 . 750 , 

7 . 000 , 

7 . 250 , 

7 . 500 , 

8.000 } 


WcMap 

= { 30 . 446 , 

30 . 533 , 

30 . 568 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

} 

NcDes= 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30.572 } 

70.000 { 







PRdes 

f 







WcMap 

= { 30 . 299 , 

30 . 413 , 

30 . 480 , 

30 . 516 , 

30 . 529 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

} 

NcDes= 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30.530 } 

80.000 { 







PRdes 

l 







WcMap 

= { 30 . 120 , 

30 . 239 , 

30 . 317 , 

30 . 368 , 

30 . 398 , 

30 . 415 , 

30 . 421 , 

30 . 421 , 

30 . 421 , 

30 . 421 , 

30 . 421 , 

30 . 421 , 

30 . 421 , 

30 . 421 , 

} 

NcDes= 

30 . 421 , 

30 . 421 , 

30 . 421 , 

30 . 421 , 

30 . 421 , 

30.421 } 

90.000 { 







PRdes 

j 







WcMap 

= { 30 . 014 , 

30 . 124 , 

30 . 201 , 

30 . 253 , 

30 . 288 , 

30 . 311 , 

30 . 325 , 

30 . 333 , 

30 . 337 , 

30 . 337 , 

30 . 337 , 

30 . 337 , 

30 . 337 , 

30 . 337 , 

} 

NcDes= . 

30 . 337 , 

30 . 337 , 

30 . 337 , 

30 . 337 , 

30 . 337 , 

30.337 } 

L 00. 000 { 







PRdes 

l 







WcMap 

= { 29 . 856 , 

29 . 948 , 

30 . 013 , 

30 . 059 , 

30 . 091 , 

30 . 113 , 

30 . 128 , 

30 . 139 , 

30 . 145 , 

30 . 149 , 

30 . 150 , 

30 . 150 , 

30 . 150 , 

30 . 150 , 

} 

NcDes= : 

30 . 150 , 

30 . 150 , 

30 . 150 , 

30 . 150 , 

30 . 150 , 

30.150 } 

L 10.000 { 







PRdes 

j 







WcMap 

= { 29 . 799 , 

29 . 870 , 

29 . 920 , 

29 . 955 , 

29 . 979 , 

29 . 997 , 

30 . 009 , 

30 . 017 , 

30 . 023 , 

30 . 026 , 

30 . 028 , 

30 . 029 , 

30 . 029 , 

30 . 029 , 

} 

30 . 029 , 

30 . 029 , 

30 . 029 , 

30 . 029 , 

30 . 029 , 

30.029 } 

} 

PgeomDesign = 2.0 { 







NcDes= 

60.000 { 







PRdes 

{ 3 . 000 , 

3 . 250 , 

3 . 500 , 

3 . 750 , 

4 . 000 , 

4 . 250 , 

4 . 500 , 


4 . 750 , 

5 . 000 , 

5 . 250 , 

5 . 500 , 

5 . 750 , 

6 . 000 , 

6 . 250 , 


6 . 500 , 

6 . 750 , 

7 . 000 , 

7 . 250 , 

7 . 500 , 

8.000 } 


WcMap 

= { 30 . 446 , 

30 . 533 , 

30 . 568 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

} 

NcDes= 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30 . 572 , 

30.572 } 

70.000 { 







PRdes 

j 







WcMap 

= { 30 . 299 , 

30 . 413 , 

30 . 480 , 

30 . 516 , 

30 . 529 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 


30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30 . 530 , 

30.530 } 
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NcDes= 80.000 { 
PRdes = *; 
WcMap = { 


} 

NcDes= 90.000 { 
PRdes = *; 
WcMap = { 


30.120, 

30.421, 

30.421, 


} 

NcDes= 100.000 { 
PRdes = *; 
WcMap = { 


30.014, 

30.333, 

30.337, 


} 

NcDes= 110.000 { 
PRdes = *; 
WcMap = { 


29.856, 

30.139, 

30.150, 


29.799, 

30.017, 

30.029, 


30. 

.239, 

30. 

.317, 

30. 

,368, 

30. 

.421, 

30. 

.421, 

30. 

,421, 

30. 

.421, 

30. 

.421, 

30. 

,421, 

30. 

.124, 

30, 

.201, 

30. 

,253, 

30. 

.337, 

30. 

.337, 

30. 

,337, 

30. 

.337, 

30. 

.337, 

30. 

,337, 

29. 

.948, 

30. 

.013, 

30. 

,059, 

30. 

.145, 

30. 

.149, 

30. 

,150, 

30. 

.150, 

30, 

.150, 

30. 

,150, 

29. 

.870, 

29. 

.920, 

29. 

,955, 

30. 

.023, 

30. 

.026, 

30. 

,028, 

30. 

.029, 

30. 

.029, 

30. 

,029, 


30. 

.398, 

30. 

.415, 

30. 

.421, 

30. 

,421, 

30. 

.421, 

30. 

.421, 

30. 

.421, 

30. 

.421 } 



30. 

.288, 

30. 

.311, 

30. 

.325, 

30. 

.337, 

30. 

.337, 

30. 

.337, 

30. 

.337, 

30. 

.337 } 



30. 

.091, 

30. 

.113, 

30. 

.128, 

30. 

.150, 

30. 

.150, 

30. 

.150, 

30. 

.150, 

30. 

.150 } 



29. 

.979, 

29. 

.997, 

30. 

.009, 

30. 

.029, 

30. 

.029, 

30. 

.029, 

30. 

.029, 

30. 

.029 } 




PgeomDesi gn . i nterp = "linear" ; 
PgeomDesi gn . extrap = "none" ; 


NcDes.interp = "linear" ; 
NcDes. extrap = "none" ; 


PRdes. interp = "linear" ; 
PRdes. extrap = "none" ; 
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compType = "Compressor"; 
compTypeFormat = "????????"; 

compAttr = { "S_map.wcMap:????.??=wcMap", "S_map.PRmap:??.???=PRmap", "S_map.effMap:??.????=effMap" , 

"S_map. NcMap :????? .???=NcMap" , "S_map. RlineMap:???.????=R/Parm" , "S_map. s_wcDes :????? .????=s_wcDes" , 
"S_map . s_PRdes : ??? . ????=s_PRdes" , "S_map . s_ef fDes : ? . ????=s_effDes" , "S_map . s_NcDes : ??? . ????=s_NcDes" 



T 3 >, LL d) 

TJ (U O) L I 

^ O a Q.4-JI— 
UCO> >. -M O 
O d) I- I- < u 
i— i— a. o. a. 3 
COPEEEO 
I--I- O O O-C 
I— 4-J U U U l/l 


</) >t LL <U 

-OdJdll-l 

^ O a Q.-MI— 
u ca >> >,+j o 
O CU I- I- < u 
I— I— Cl Cl Q. 2 
CA -M E E E O 
|--i- O O O-C 
i— 4 -> U U U l/> 


l/l >, LL d) 
■a di di s- i 
^ O CL Q.4-JI— 
U CA >, >. -M o 
O d) I- I- < u 
I— I— Cl Cl Cl 3 
Cfl +J E E E O 
|--i- O O O-C 
i— 4 -> U U U i/> 


"D dJ d) C I 

^ O a Q.-MI— 
U Cfl >> >,-M o 
o <u i— I— < u 
I— I— Cl Cl Cl 2 
Cfl -M E E E O 
I-t O O O-C 
i— 4-> U U U l/l 


l/l -M 

■- S-4- ■- 
: dJ rd = 

i/i 4-1 x: i/i 
■M 4-J l/l ~0 
U-i-: cd 

3 1 — O l/l 

1 -~0 Q- II i — 4- 1 
r i/l = C 


I— E 4-> E a. 
C O E O O E 
•r 4-J O-Q -M O 
: 4-14-1 • 4-< U 

O 4-< i/> O 
II XI O S-X14- 
•XI <U • O 
S— l/l • 4-* l/l 
O 4-> l/l 4-/ 4-* T5 
SZ dJ 4-* -1— M— C 
Ur- Ur- n 3 d) 
c c 3 clx: \ 
-V 03 -I- T3 l/l l/l\ 


+ 3 : : 3 N 

+ 

s- r- xi r- 
+ aioc'-- 3 

+ (JIN || 4 — ’ 

IA 4 JNM- 

--'-Y-' S_ 3NM 

(U OC'-C'-- 
■■Ifl+J'DNN - 
- - "a c cu r- ■ : 

-• d> -i— cu c'-- c'-- c'-- 

-- CU r- C'-C'-C'-- 
+ i— I Cfl 0 s " •■C'" 
+ Cfl (UN 3 N 
\ di i/i air- "a ■ 
\ aid cd ('-• i— ('-■ 
cd lu -m c'-- cfl C'-- 

4 -> LU l/lr UN 
Vl/I J L fd C'-- 

s_ cfl cu || s_ ■■ 
o. <u c -m 4-3 
3 4 CP: : 

O CZ H fd 


■a jlC ll o ll 

d) U > LL 
d> o "a d) cu l. 
r— I— O a. Q. -M 

Cfl Cfl Cfl >> >. -M 

I— CU I— I— < 

i— Q. Q. CL 
U O -M E E E 

O U-r- O O O 

i— 4 -> U U U 
Cfl C 


4-1 fv.. c^. (n.. 

i/i i s_ cv.. ■ r- 

-a oc'- c^- c'-- 

d) l/l Q- C'" C'" ■ 

dJ Q (JN- •■C'-- 
i — LU 3 N 3 N 

cfl lu o= u r- 

4-J _IT3 fd C'-- 
d) Cfl d) II L ■ ■ 
I- C di 4-3 

3 : Cfl fd 


U > LL 
O U CU CU L. 
I— O CL Q.-M 
Cfl Cfl ^ >. -M 

I— <u |— I— < 

I— I— Q. Q. Q. 

O 4 -< E E E 

U ■!— o o o 

4 -< U U U 


l/l l/l L. U 

-a —i u cd 

CU CO a. 

CU LL CU l/l 


an- o 'd- 

fd (U 4 - i— II 

■M -M O Cfl || 

LO 4-» 4-* 4-* 

i— 3 T 3 xrr 

diu c cu 4 -> ai 


NASA/TM— 2007-2 14690 


32 



anchor = "burners"; 
burners . bottom = "mixers" 
mixers . bottom = "nozzles' 
//end of Cblock 
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TURBOMACHINERY MAP 
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Appendix C 
Turbomachinery Maps 

In the performance analysis of gas turbine engines, it is very desirable to have relationships that 
describe the variation in mass flow rate, rpm, pressure ratio, and temperature ratio for each turbomachine 
in the model. These characteristics are often calculated by complex codes using velocity diagrams, fluid 
mechanics, empirical correlations, and even computational fluid dynamics. The physics of 
turbomachinery flow simplify when corrected (or equivalent) airflow is used rather than absolute flow 
and corrected ( or equivalent) speed is used for absolute lpm. The calculation of the relationship between 
corrected flow, corrected speed, pressure ratio, and temperature ratio (or efficiency) for a turbomachine 
results in a collection of data referred to as a “map”. 

Figure 2 shows a compressor map, normally plotted as pressure ratio (across the machine) versus 
corrected airflow into the machine. For multiple values of corrected rpm the relationship between 
pressure ratio and flow is shown and the line that connects the upper left of each speed characteristic is 
called the surge line. Each speed characteristic also has a lower limit: the line becomes vertical as the 
machine chokes and no further mass flow can be passed. Finally, values of constant adiabatic efficiency 
are plotted; these contours resemble, roughly speaking, a series of concentric ovals and are often referred 
to as “efficiency islands”. 

Figure 3 shows the same map with some of the data values identified. The application of such a 
specific compressor map to a complete engine system is quite limited; however, if the map could be 
scaled its usefulness would be dramatically increased. NPSS allows the map to be scaled up or down in 
pressure ratio, efficiency, flow, or any combination of the three. Assume that the design point is identified 
as 100 percent speed and 25 percent surge margin, resulting in an unsealed pressure ratio of 2.0, corrected 
airflow of 500 lb m /s, and adiabatic efficiency of 87 percent. Further assume that the desired values of 
pressure ratio, flow, and efficiency at that point are 1.80, 475 lb m /s, and 87.5 percent respectively. The 
ratio of the desired value to the unsealed value of each parameter determines the scale factor applied to 
the map (one scale factor for each parameter). 
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Figure 4. — Scaled compressor map. 

Figure 4 is the scaled compressor map. The scale factors are calculated during the design point 
analysis and used during the performance analysis from the following equations: 
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flow scale factor = WCdesired = 475 lb/s = 0,950 
C unsealed 500 lb/ S 

PR scale factor = ^desired" 1 ^ 18Q ~ 1 ^ 0 800 
PR unsealed — 1 2.00 — 1 

efficiency scalefactor= P desired _ 0-875 
T unsealed 0.870 

Note that the limit for pressure ratio is 1, not 0, because the compressor/turbine requires/performs no 
work at that pressure ratio; this is why 1 is subtracted from both values in the pressure ratio scale factor 
equation. These scale factors are applied to every point in the map: thus all values of corrected flow are 
multiplied by 0.95, all efficiencies are increased by 0.6 percent, and pressure ratio minus one is reduced 
20 percent as shown by the two points in figures 3 and 4. 

There are two peculiarities of the compressor corrected speed characteristic because of its shape 
(resembling a candy cane): first, it is possible to have a particular value of pressure ratio and corrected 
speed occur at two values of corrected flow; second, pressure ratio is indeterminate when a particular 
value of corrected speed and flow occur on the vertical part of the speed characteristic. Because of this, 
compressor maps require an second independent parameter (the first is corrected speed) to uniquely 
determine the flow, PR and efficiency: for one commonly used format of compressor map in NPSS that 
parameter is called the R-line. 

Figure 5 shows the compressor map with R-lines included. Typically the first R-line is given a value 
of 1.0 and coincides with the surge line. The remaining R-lines are roughly parallel to the surge line with 
increasing R-line value corresponding to higher surge margin. These R-lines have no physical meaning; 
they are simply mathematical constructs orthogonal to corrected speed so that any value of R-line and 
corrected speed will uniquely position the compressor operating point. In other words, given R-line value 
and corrected speed, the compressor pressure ratio, corrected flow, and adiabatic efficiency are known, 
and further, the derivatives of these parameters can be calculated with respect to R-line and corrected 
speed which is critical for computational convergence. 

With regard to efficiency and corrected flow, the scaling of a turbine map (shown in fig. 6) is the 
same as with compressor maps. Since pressure ratio and speed uniquely determine corrected flow and 
efficiency, no R-lines are required for turbine maps. The scaling of pressure ratio, however, is slightly 
different: 


turbine PR scale factor = PRdesired 1 = 4 ' 0Q 1 =1.163 

PRscaled - 1 3.58 — 1 

In this example the desired turbine pressure ratio (the unsealed pressure ratio used to look up the turbine 
corrected flow and efficiency) is 4.00 while the actual turbine pressure ratio is 3.58, resulting in a scale 
factor of 1 . 1 63. NPSS will scale maps by any factor; it is up to the user to determine the validity of the 
resulting map. Certainly a map representing a single-stage axial fan should not be scaled to a design 
pressure ratio of 10. 

Lastly, some maps have an additional input variable such as inlet guide vane (IGV) setting or nozzle 
vane angle. In this case there are multiple plots similar to figure 2 or figure 6 that comprise the 
turbomachine map, one such plot for each value of the additional variable. These are referred to as 
“stacked maps” or “3-D maps” and the third-dimensional variable is named “alpha” for compressors and 
‘parmGeom’ for turbines. 
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